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ABSTRACT 

The structure and properties of the diffuse interstellar medium (ISM) on small scales, 
sub-au to 1 pc, are poorly understood. We compare interstellar absorption-lines, ob¬ 
served towards a selection of O- and B-type stars at two or more epochs, to search 
for variations over time caused by the transverse motion of each star combined with 
changes in the structure in the foreground ISM. Two sets of data were used: 83 VLT- 
UVES spectra with approximately 6 yr between epochs and 21 McDonald observatory 
2.7m telescope echelle spectra with 6 - 20 yr between epochs, over a range of scales 
from ^ 0 - 360 au. The interstellar absorption-lines observed at the two epochs were 
subtracted and searched for any residuals due to changes in the foreground ISM. Of 
the 104 sightlines investigated with typically hve or more components in Nai D, possi¬ 
ble temporal variation was identified in five UVES spectra (six components), in Ca ii, 
Ca I and/or Nai absorption-lines. The variations detected range from 7% to a fac¬ 
tor of 3.6 in column density. No variation was found in any other interstellar species. 
Most sightlines show no variation, with 3cr upper limits to changes of the order 0.1 - 
0.3 dex in Can and Nai. These variations observed imply that fine-scale structure is 
present in the ISM, but at the resolution available in this study, is not very common 
at visible wavelengths. A determination of the electron densities and lower limits to 
the total number density of a sample of the sightlines implies that there is no striking 
difference between these parameters in sightlines with, and sightlines without, varying 
components. 

Key words: stars: early-type - ISM: abundances - ISM: clouds - ISM: general - 
ISM: structure - 


1 INTRODUCTION 

The properties of the interstellar medium (ISM) on small 
scales, from sub-astronomical unit (au) to 1 parsec, remain 
poorly understood. The first indication that au-scale struc¬ 
tures were present in H l was provided by |Diamond et al.| 
119891, who used Very Long Baseline Interferometry (VLBI) 
radio observations to identify changes in the 21 cm absorp¬ 
tion towards three extragalactic radio sources as a function 


of resolution. This suggested foreground structures with lin¬ 
ear dimensions of only ~25 au. Further absorption-line VLBI 
measurements in H l (Faison et al.|1998 Brogan et al.|2005 


Lazio et al.|2009] |Roy et al.|2012[) and molecules ( |Marscher,| 

Moore fc Bania||1993 Liszt &: Lucas|[^000l also found tiny- 


scale structures, although with relatively small plane-of-sky 
covering fractions of ~10 per cent. Similar results have been 
seen via multi-epoch Hi absorption-line measurements to¬ 


wards pulsars (e.g. Frail et al. 1994 Johnston et al. 2003 


|Stanimirovic et al.|2010[) and indicate that, in Hi, tiny-scale 
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variations are a common phenomenon but with a low prob¬ 
ability of intercepting any given line of sight. 

Optical and UV observations have also been used to 
detect temporal and spatial variation in interstellar absorp¬ 
tion. Most have relied upon multi-epoch ultra-high resolu¬ 
tion optical observations in Nal and Call (e.g [Hobbs et a~ 


1991[ [Blades et al.| |1997l |Price, Crawford fc Barlow||2000j 


Meyer et al.|2012|, although some authors have also de- 


tected variations in Kl, Cal and/or Fel 1 Crawford et al. 

2000 Smith et al. 2013 Galazutdinov et al.[|2013l. Binary 

systems (e.g. Meyer 1990 Cordiner et al. 

|2013|) and star 

clusters (e.g. Andrews, Meyer & Lauroesch 

[20011 van Loon 


|et al.|200^ have also been used to study the fine-scale struc¬ 


ture of the ISM. Watson & Meyer (1996 \ found differences 


in Na I towards the individual members of all 17 of the 
binary or multiple star systems they observed, with separa¬ 
tions 480-29000 au. These observations imply a population 
of dense au-scale clouds in the ISM or changes in ionization 
or depletion over time, but are restricted to trace neutral 
species. Boisse et al. (20131 studied the spatial distribution 


of CH"*", CH and CN on scales ranging 1-20 au for a single 
sightline towards Zeta Per, finding the neutral species to be 
uniformly distributed and a CH^ excess implying its pro¬ 
duction in very small localized active regions of only a few 
au. Data for UV absorption-lines can probe dominant ions 
and provide further information on the physical conditions; 
time-variable UV lines have been found by e.g. |Danks et al.| 
| |2001[ ) and |Welty| ( |2007[ ) . Further examples of time-variable 
interstellar lines are cited in the reviews by Crawford (20031 


and Lauroesch (20071. 


We note it is difficult to decipher between apparent spa¬ 
tial and temporal variations in the foreground ISM. The 
main questions concerning these variations are as follows. 
What are the physical conditions in these structures? What 
is their formation mechanism? How stable are they and 
how do the structures seen in the optical relate to the Hi 
absorption-line measurements? 

This is one of a series of papers which use data on early- 


Bagnulo et al. 

20031 

jto study the ISM. In 

Hunter et al. 


(2006 Paper I) the near-UV interstellar lines of Nal, Till 
and Call were observed towards 74 O- and B-type stars, 
and basic correlations in the observed column densities were 
investigated. [Smoker et al.| | |2007[ Paper H) combined these 
data with H l 21 cm observations to estimate distances to 
intermediate- and high-velocity clouds, while [Smoker et al.| 
| [2011[ Paper HI) used POP observations of stars in open 
clusters (IC2391, NGC6475 and Praesepe) combined with 
archival multi-epoch data to search for small- (pc) and tiny- 
(au) scale structure of the ISM. They searched for variations 
in Till, Call K, Nal D and Kl absorption; a total of three 
possible examples of tiny-scale structure were found. 

Here we present UVES observations at additional 
epochs of the same POP O- and B-type lines of sight, cou¬ 
pled with additional data taken at McDonald observatory to 
search for time-variable absorption. The time between the 
two epochs was ~6 yr for the UVES POP sample and ~ 6 
- 20 yr for the McDonald data. Due to the proper motion 
of the stars, it is hence possible to study changes in their 


^ http://www.eso.org/uvespop 


lines of sight on scales of ~ 0 - 360 au. We aim to extend 
the sample size of twin epoch observations to determine just 
how ubiquitous the small scale structure of the ISM is. To 
our knowledge this is the largest twin-epoch spectroscopic 
survey to date. 

This paper is organized as follows. Section describes 
the archival data used, the new observations, and the data 
reduction, while Section [^contains the main results. In Sec¬ 
tion]^ we discuss the upper limits and possible time-varying 
interstellar line strength in terms of different models of the 
ISM. Finally, Section presents a summary of our results 
and suggests avenues for future research. 


2 OBSERVATIONS, ARCHIVAL DATA, AND 

DATA REDUCTION 

2.1 UVES data 

Our study uses a combination of archival and new observa¬ 
tions. Archival data from the ESO POP (ESO programme 
ID 266.D-5655(A), [Bagnulo et al.|200^ were used to provide 
a catalogue of high resolution {R= ~ 80, 000), covering 

the wavelength range 304 - 1040 nm for about 400 stars ob¬ 
tained with the upgraded UVES spectrometer ( [Dekker et al.[ 
2000 Smoker et al.|200^ of the ESO Very Large Telescope. 
The signal-to-noise ratio (S/N) of the stellar spectra is high, 
ranging 200 - 800, and averaging 520 around Nal, and rang¬ 
ing 100 - 800, averaging 330 around Ca ll. The POP sur¬ 
vey contains a total of 98 O- and B-type stars, suitable for 
use as background sources to study the ISM. These sight¬ 
lines form the first epoch of data for most stars in our sam¬ 
ple, and were observed between 2001 and 2003. The merged 
versions of the spectra were taken from the POP website 
((http://www.eso.org/uvespop Bagnulo et al.[[^003 l). 

New second-epoch observations of the same O- and B- 
type stars were taken during 2008 (ESO programme ID 
081.C-0475(A)) on the same UVES instrument to allow di¬ 
rect comparison between the two epochs at the same spectral 
resolution. It is important to note that the measured reso¬ 
lution of the UVES data can chanre by up to 10 % over 
several yr, depending on the season R We have investigated 
the effects of these small resolution changes on the column 
densities calculated and found that in the majority of cases. 


the differences are within the error estimates (see Sect. 2.51. 
The (390-1-580) setting was used giving wavelength coverage 
from 328-450 nm and 475-681 nm, with a gap from 578-581 
nm where the two CCDs in the red mosaic are joined. For 
some sightlines where changes were detected, we obtained 
third epoch data in 2013 using the same instrumental setup 
(ESO programme ID 092.C-0173(A)). All of these data were 
reduced using the ESO UVES pipeline ( [Ballester et al.|2000[ ) 
and the multiple exposures at each epoch co-added using 
IRAI0 to bring the epoch 2 and 3 data to the same point of 
reduction as epoch 1. As the epoch 1 data were taken from 


^ http://www.eso.org/observing/dfo/quality/UVES/reports/ 
HEALTH/trend_report_ECH_RESQLUTIDN_DHC_HC. html 

IKAh' is distributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association of Universities 
for Research in Astronomy (auRA) under cooperative agreement 
with the National Science Foundation. 
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the POP website, they had been reduced using an earlier 
version of the pipeline than the epoch 2 and 3 data. To in¬ 
vestigate the effects of the different pipelines on the spectral 
resolution, a subset of epoch 1 data (including some show¬ 
ing possible variation) were re-reduced using a more recent 
pipeline. No significant resolution changes were discovered 
and so the merged version of the data from the POP website 
was deemed sufficient for comparison. 

From this stage onwards all sets of UVES data were 
treated in exactly the same manner. Each spectrum was 
analysed in idl. Interstellar absorption-lines of interest 
(listed in Table were isolated and each line from epoch 1 
was cross correlated with its equivalent from epoch 2 or 3 to 
correct for any velocity shifts between the ISM lines in the 
two epochs, and to bring both epochs into kinematic local- 
standrad-of-rest frame. The two epochs were then overlaid 
so that any obvious changes in the line profiles could be 
identihed by eye. Spectra were renormalized to remove the 
stellar continua and the Nal D spectra corrected for telluric 
absorption features where necessary, by dividing by a scaled 
telluric spectrum created from our set of observations with 
regions free of Na absorption or from Cerro Paran al Sky 
Model SKYCALcP] ( |N o11 et al.|20T2 Jones et al.|2013 1. In all 
cases where possible variations were identified,they are seen 
in both the Na Di and Na D 2 lines, which are affected differ¬ 
ently by telluric contamination, and so inconsistent correc¬ 
tion of telluric lines could not account for the variation seen, 
epoch 1 spectra were subtracted from epoch 2 or 3 for each 
absorption feature to search for any residuals that would 
imply variations in the line of sight between epochs. It was 
also necessary to compare these residual plots to residuals 
caused solely by the slight differences in resolution in the 
two epochs of data due to small, but inevitable changes in 
the instrumental set-up between epochs-exposures. 


2.2 Welty et al. data 


In some cases, the earlier epoch 1 data used for comparison 


were taken from Welty, Hobbs & Kulkarni (1994 Henceforth 


refered to as WHK) and Welty, Morton & Hobbs 11996 


Henceforth WMH). WHK observations of Na l Di only, were 
made in 1987-89 on a coude spectrograph of the 2.7m tele¬ 
scope at McDonald observatory. WMH observations of Ca ll 
were obtained between 1993 and 1995 with the 0.9m coude 
feed telescope and 2.1m coude spectrograph of the Kitt Peak 
National Observatory, and with the Ultra-High Resolution 
Facility (UHRF) on the 3.9 m Anglo-Australian Telescope 
in 1994. All of these epoch 1 data have very high S/N rang¬ 
ing 100-700 averaging around 300, and spectral resolution 

of R ~ 214,000-938,000. Twenty-one sightlines had two 

epochs available so the McDonald data set is therefore signif¬ 
icantly smaller than that from UVES, but has nearly twice 
the spectral resolution and a significantly longer interval be¬ 
tween epochs. The WHK Nal spectra have higher resolution 
and so were degraded using a Gaussian of appropriate reso¬ 
lution, for comparison with our later observations. Data were 
reduced using standard methods within iraf. Spectra were 
renormalized to remove stellar continua and telluric correc- 


^ http://www.eso.org/observing/etc/ 
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Figure 1. Histogram of transverse distance moved for the stars 
in the UVES and McDonald samples between the two epochs. 


tions were applied where necessary in the same manner as 
described above for the UVES data set. 


2.3 McDonald Data 


Additional observations were taken at the W.J. McDonald 
Observatory with the Tull coude echelle spectrograph on 
the 2.7m Harlan J. Smith telescope ( Tull et al.|1995 i during 
May-June 2008 and July 2014. Spectra were acquired at 
two different wavelength settings. One setting provided the 
principal interstellar lines Ca ll K and Ca l at 4226A, as well 
as the potential for detecting lines from K l at 4044A and 
molecules CH, CH’^ and CN. A second setting covered the 
interstellar Nal D lines. Both settings provided a resolving 
power R ~ 140, 000 — 160, 000. 


The hnal stellar sample is outlined in Table 2. The 
galactic coordinates, V-band magnitude, B-V colour and 
E(B-V) reddening (Hunter et al. 2006), proper motions, par¬ 
allaxes (van Leeuwen 2007) and spectral types (taken from 
the POP web site (Bagnulo et al. 2003), WHK and WMH) 
are listed for each sightline. 


2.4 Stellar distances and proper motions 


Stellar distances, (B — V) and E{B — V) values were, in 
most cases, taken from [van Leeuwen| (|2007[) and |Hunter| 


et al. (20061, taking the most accurate measurement in each 


case. For some sightlines neither Hunter et al. (20061 nor 
[van Leeuwen| ( |2007[ ) gave distance estimates and so values 
were taken from the literature. These cases are referenced 
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Table 2. The stellar sample. Columns are as follows: the HD number of the star; the galactic coordinates I and b; the V band magnitude, 
B — V colour and reddening E(B — V") | |Hunter et al.|2006[ |; the proper motions in right ascension (pa) and declination (p^), the parallax 
(tt) (from j^n Leeuwen] l |2007| ) the spectral types and the instrument used to obtain our data set. 


Star 


1 

(deg) 

b 

(deg) 

V 

(mag) 

B-V 

(mag) 

E(B - V) 
(mag) 

pa 

(mas yr“^) 

P5 

(mas yr~^) 

TT 

(mas) 

Spectral 

type 

Data set 

HD 

10840 

291.7609 

-54.9427 

6.79 

-0.08 

- 

26.93 ± 0.29 

-5.47 ± 0.34 

5.55 ± 0.36 

B9p 

UVES 

HD 

29138 

297.987 

-30.5423 

7.19 

-0.09 

0.1 

-2.53 ±0.36 

8.93 ±0.45 

0.16 ± 0.34 

Bllab 

UVES 

HD 

30677 

190.1808 

-22.2169 

6.84 

-0.02 

0.17 

-3.13 ±0.57 

-4.00 ±0.39 

0.89 ± 0.53 

BlH/HI((n) 

UVES 

HD 

33328 

209.1402 

-26.6856 

4.25 

-0.18 

0.01 

0.25 ±0.16 

-1.97 ±0.11 

4.02 ± 0.18 

B2IVne 

UVES 

HD 

43285 

203.4236 

-05.1276 

6.05 

-0.11 

0.03 

-0.40 ±0.42 

-17.42 ±0.32 

4.99 ± 0.43 

B6Ve 

UVES 

HD 

45725 

216.6614 

-08.2139 

4.6 

-0.1 

0.08 

-7.0 ±1.2 

-5.0 ±0.9 

- 

B3Ve 

UVES 

HD 

46185 

221.9706 

-10.0784 

6.79 

-0.16 

-0.02 

-4.18 ±0.57 

-0.31 ±0.48 

2.63 ± 0.55 

B2/B3II: 

UVES 

HD 

46462 

245.4745 

-19.6548 

7.56 

-0.10 

0.52 

-8.81 ±0.41 

9.43 ±0.57 

2.93 ± 0.52 

B9 

UVES 

HD 

47116 

233.174 

-13.9282 

7.71 

0.01 

0.52 

0.06 ±0.44 

5.43 ±0.49 

1.53 ± 0.61 

B9p 

UVES 

HD 

48099 

206.2096 

+00.7982 

6.37 

-0.08 

0.21 

0.84 ±0.41 

2.5 ±0.31 

1.17 ± 0.41 

07V 

UVES 

HD 

49131 

240.5001 

-14.7267 

5.8 

-0.19 

0 

-2.49 ±0.41 

3.80 ±0.53 

1.72 ± 0.54 

B2HI 

UVES 

HD 

49333 

231.3595 

-10.3419 

6.08 

-0.19 

-0.15 

-14.53 ±0.43 

5.12 ±0.51 

4.14 ± 0.51 

B7p 

UVES 

HD 

50896 

234.7568 

-10.0832 

6.74 

-0.06 

0.17 

-3.76 ±0.33 

4.35 ±0.56 

0.72 ± 0.62 

WN5 

UVES 

HD 

51876 

228.0469 

-05.7456 

7.13 

-0.08 

-0.07 

-0.27 ±0.48 

8.70 ±0.50 

3 ± 0.59 

BOH 

UVES 

HD 

52918 

218.0123 

+00.6139 

4.99 

-0.2 

0.03 

-5.05 ±0.26 

2.24 ±0.19 

2.68 ± 0.22 

BIV 

UVES 

HD 

58343 

231.0872 

-00.2141 

5.2 

-0.04 

0.17 

-4.95 ±0.47 

-20.21 ±0.45 

3.48 ± 0.58 

B2Vne 

UVES 

HD 

58377 

242.1965 

-06.2584 

6.81 

-0.17 

-0.02 

-7.14 ±0.32 

5.69 ±0.48 

1.42 ± 0.47 

B5IV 

UVES 

HD 

58978 

237.4106 

-02.9978 

5.61 

-0.13 

0.06 

-7.76 ±0.25 

4.82 ±0.34 

2.61 ± 0.33 

Bill 

UVES 

HD 

60498 

247.1984 

-06.6420 

7.35 

-0.12 

0.03 

-8.63 ±0.29 

4.18 ±0.34 

4.19 ± 0.34 

B4HI 

UVES 

HD 

61429 

240.6523 

-01.8360 

4.7 

-0.11 

-0.01 

-3.03 ±0.13 

-7.42 ±0.22 

5.27 ± 0.26 

B8IV 

UVES 

HD 

62714 

268.5841 

-15.6252 

7.25 

-0.08 

0.05 

-14.30 ±0.51 

7.15 ±0.49 

4.1 ± 0.41 

B7Vp 

UVES 

HD 

64972 

245.0907 

-00.0164 

7.18 

-0.11 

-0.05 

-9.1 ±1.5 

6.8 ±1.3 

- 

B8p 

UVES 

HD 

67536 

276.1126 

-16.1399 

6.24 

-0.09 

0.12 

-5.36 ±0.27 

11.43 ±0.21 

2.72 ±0.22 

B2.5Vn 

UVES 

HD 

68761 

254.3713 

-01.6188 

6.53 

-0.09 

0.13 

-4.44 ±0.26 

5.71 ±0.29 

1.63 ± 0.34 

B0.5HI 

UVES 

HD 

72067 

262.0769 

-03.0778 

5.83 

-0.16 

0.05 

-7.00 ± 0.59 

6.98 ± 0.41 

2.28 ± 0.49 

B2Vne 

UVES 

HD 

74966 

258.0811 

+03.9290 

7.43 

-0.14 

0.01 

-10.03 ± 0.36 

6.98 ± 0.38 

2.24 ± 0.52 

B4IV 

UVES 

HD 

76131 

273.5774 

-07.2675 

6.69 

-0.05 

0.08 

-16.52 ± 0.30 

6.16 ± 0.25 

2.38 ± 0.29 

B6HI 

UVES 

HD 

76341 

263.5348 

+01.5210 

7.17 

0.22 

0.49 

-5.16 ± 0.40 

3.49 ± 0.44 

0.75 ± 0.48 

09.7Iab 

UVES 

HD 

83312 

284.7656 

-10.4703 

7.47 

-0.01 

0.15 

-13.47 ± 0.48 

7.24 ± 0.45 

1.93 ± 0.4 

B4;psh 

UVES 

HD 

88661 

283.0814 

-01.4804 

5.76 

-0.12 

0.08 

-10.56 ± 0.37 

6.81 ± 0.35 

2.56 ±0.37 

B2IVe 

UVES 

HD 

89587 

279.8352 

+05.1930 

6.87 

-0.13 

0.03 

-10.99 ± 0.33 

6.46 ± 0.31 

0.89 ± 0.38 

B3HI 

UVES 

HD 

90177 

285.1466 

-01.9838 

7.57 

0.86 

1.07 

-6.87 ± 0.99 

2.33 ± 0.77 

1.69 ± 0.82 

B2ev 

UVES 

HD 

90264 

289.0967 

-08.1064 

4.95 

-0.12 

-0.01 

-22.39 ± 0.20 

11.48 ± 0.18 

8.12 ± 0.18 

B8p 

UVES 

HD 

90882 

248.4149 

+44.5898 

5.18 

-0.04 

0.01 

-48.86 ± 0.22 

-13.43 ± 0.19 

10.13 ± 0.23 

B9.5V 

UVES 

HD 

92740 

287.1742 

-00.8475 

6.42 

0.08 

0.08 

-7.71 ± 0.36 

2.86 ± 0.3 

0.08 ± 0.34 

WN7A 

UVES 

HD 

94910 

289.1829 

-00.6953 

7.09 

0.49 

0.08 

-5.89 ± 0.83 

2.45 ± 0.76 

0.19 ± 0.7 

B2;pe 

UVES 

HD 

94963 

289.7649 

-01.8088 

7.14 

-0.1 

0.19 

-7.9 ± 1.3 

2.8 ± 1.3 

- 

06.5HIe 

UVES 

HD 

100841 

294.4727 

-01.3974 

3.12 

-0.04 

0.03 

-33.41 ± 0.33 

-7.08 ± 0.29 

7.77 ±0.34 

B9HI 

UVES 

HD 

105056 

298.9454 

-07.0551 

7.42 

0.09 

0.26 

-4.34 ± 0.44 

-1.28 ± 0.4 

0.21 ± 0.49 

ON9.7Ia 

UVES 

HD 

105071 

298.2376 

-03.0911 

6.31 

0.13 

0.19 

-5.76 ± 0.43 

-0.77 ± 0.32 

0.05 ± 0.41 

B6Ia/Ib 

UVES 

HD 

106068 

298.5069 

-00.4152 

5.95 

0.23 

0.26 

-5.07 ± 0.30 

-0.52 ± 0.3 

0.08 ± 0.34 

B8Ia/Iab 

UVES 

HD 

109867 

301.7113 

-04.3511 

6.26 

0 

0.19 

-4.54 ± 0.31 

-1.09 ± 0.28 

-0.54 ± 0.32 

Blla 

UVES 

HD 

112272 

303.4864 

-01.4947 

7.39 

0.64 

0.84 

-8.99 ± 0.42 

-8.10 ± 0.34 

1.4 ± 0.53 

BO.SIa 

UVES 

HD 

112842 

304.0558 

+02.4766 

7.08 

0.15 

0.25 

-4.24 ± 0.45 

-0.95 ± 0.38 

0.31± 0.52 

B4IV 

UVES 

HD 

113904 

304.6745 

-02.4907 

5.69 

-0.11 

0.11 

-4.26 ± 0.39 

-2.18 ± 0.48 

0.26 ± 0.48 

WC5/B0Ia 

UVES 

HD 

115363 

305.8832 

-00.9683 

7.82 

0.5 

0.69 

-8.48 ± 0.52 

-0.67 ± 0.54 

0.87 ± 0.68 

Blla 

UVES 

HD 

115842 

307.0804 

+06.8343 

6.04 

0.21 

0.41 

-3.04 ± 0.31 

3.93 ± 0.36 

0.65 ± 0.44 

B0.5Ia 

UVES 

HD 

120709 

317.2822 

+28.1893 

4.53 

-0.14 

0.01 

-34.54 ± 0.78 

-28.21 ± 0.62 

9.49 ± 0.89 

B5Ip 

UVES 

HD 

123515 

315.1159 

+09.5089 

5.96 

-0.04 

0.03 

-32.25 ± 0.32 

-11.63 ± 0.26 

3.98 ± 0.39 

B8p 

UVES 

HD 

125823 

321.5656 

+20.0226 

4.40 

-0.17 

-0.05 

-24.15 ± 0.13 

-21.90 ± 0.13 

7.13 ± 0.16 

B8p 

UVES 

HD 

133518 

323.0822 

+05.4606 

6.40 

-0.11 

0.08 

-2.83 ± 0.40 

-12.52 ± 0.48 

2.21 ± 0.44 

B2IV 

UVES 

HD 

136239 

321.2281 

-01.7448 

7.87 

0.77 

0.94 

-5.70 ± 0.99 

-2.20 ± 0.9 

0.09 ± 0.9 

B5p 

UVES 

HD 

137509 

315.2529 

-12.1253 

6.89 

-0.13 

-0.02 

-16.08 ± 0.28 

-14.64 ± 0.4 

5.11 ± 0.38 

B8p 

UVES 

HD 

137753 

325.9312 

+03.3436 

6.7 

-0.03 

0.09 

-7.79 ± 0.39 

-8.97 ± 0.38 

3.3 ± 0.42 

B7IV 

UVES 

HD 

142301 

347.1237 

+21.5122 

5.86 

-0.07 

0.04 

-13.16 ± 0.60 

-25.55 ± 0.46 

6.33 ± 0.43 

B8p 

UVES 
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Table 2. continued 


Star 


1 

(deg) 

b 

(deg) 

V 

(mag) 

B-V 

(mag) 

E(B-V) 

(mag) 

(mas yr“^) 

^5 

(mas yr~^) 

TT 

(mas) 

Spectral 

type 

Data set 

HD 

142758 

325.3065 

-04.2765 

7.08 

0.14 

0.31 

-5.38 ± 0.48 

-5.33 ± 0.47 

-0.76 ± 0.51 

Bl.SIa 

UVES 

HD 

142983 

356.3867 

-1-28.6317 

2.29 

-0.09 

0.14 

-12.44 ± 0.25 

-16.73 ± 0.21 

6.97 ± 0.24 

B8Ia/Iab 

UVES 

HD 

143448 

324.5485 

-05.9714 

7.30 

-0.00 

0.16 

-3.88 ± 0.54 

-2.72 ± 0.55 

1.54 ± 0.58 

B3IVe 

UVES 

HD 

145482 

348.1166 

-1-16.8354 

4.57 

-0.15 

0.06 

-10.38 ± 0.18 

-23.94 ± 0.14 

6.81 ± 0.16 

B2V 

UVES 

HD 

145792 

351.0096 

-1-19.0290 

6.42 

0.08 

0.14 

-8.89 ± 0.71 

-20.23 ± 0.47 

6.96 ± 0.65 

B6p 

UVES 

HD 

145842 

334.7002 

-1-02.5325 

5.12 

-0.10 

0.01 

-34.39 ± 0.25 

-44.88 ± 0.22 

8.46 ± 0.26 

B8V 

UVES 

HD 

148184 

357.9328 

-1-20.6766 

4.28 

0.31 

0.52 

-5.41 ± 0.31 

-21.12 ± 0.24 

6.21 ± 0.23 

B2Vne 

UVES 

HD 

148379 

337.2456 

-1-01.5757 

5.36 

0.44 

0.61 

-4.54 ± 1.06 

-1.51 ± 1.27 

1.79 ± 0.7 

B1.5Iae 

UVES 

HD 

148688 

340.7197 

-1-04.3474 

5.33 

0.27 

0.46 

-4.35 ± 0.37 

-1.37 ± 0.29 

1.2 ± 0.33 

Blla 

UVES 

HD 

148937 

336.3661 

-00.2181 

6.77 

0.24 

0.09 

1.51 ± 0.86 

-3.57 ± 0.77 

2.35 ± 0.79 

06.5e 

UVES 

HD 

152003 

343.3324 

-1-01.4077 

6.99 

0.29 

0.53 

-2.5 ± 1 

-2.5 ± 1 

- 

09Iab 

UVES 

HD 

152235 

343.3111 

-1-01.1041 

6.34 

0.42 

0.61 

-1.07 ± 0.56 

-3.30 ± 0.39 

1.05 ± 0.57 

B0.7Ia 

UVES 

HD 

154811 

341.0616 

-04.2191 

6.93 

0.4 

0.46 

-1.13 ± 0.62 

-3.86 ± 0.43 

2.34 ± 0.69 

OC9.7Iab 

UVES 

HD 

154873 

341.3454 

-04.1096 

6.70 

0.28 

0.47 

1.47 ± 0.72 

-4.83 ± 0.57 

2.16 ± 0.72 

Bllb 

UVES 

HD 

155416 

348.7532 

-1-00.6515 

6.90 

0.3 

0.34 

-0.26 ± 0.59 

-2.32 ± 0.27 

0.9 ± 0.54 

B7II 

UVES 

HD 

155806 

352.5858 

-1-02.8683 

5.61 

-0.06 

0.23 

0.92 ± 0.36 

-3.05 ± 0.19 

0.89 ± 0.4 

07.5Ve 

UVES 

HD 

156385 

343.1641 

-04.7627 

6.92 

0.05 

0.35 

1.75 ± 0.61 

-3.89 ± 0.3 

-0.26 ± 0.57 

WC7 

UVES 

HD 

157038 

349.9548 

-00.7940 

6.41 

0.64 

0.81 

-1.33 ± 0.52 

-0.07 ± 0.25 

0.51 ± 0.42 

Bl/B2Ia 

UVES 

HD 

163745 

350.5576 

-08.7882 

6.5 

-0.09 

-0.01 

1.43 ± 0.52 

-3.77 ± 0.21 

2.63 ± 0.4 

B5II 

UVES 

HD 

163758 

355.3605 

-06.1042 

7.32 

0.02 

0.34 

3.3 ± 1.3 

-4.8 ± 1.2 

- 

06.5Iaf 

UVES 

HD 

163800 

7.052 

-1-00.6878 

7.01 

0.31 

0.6 

-0.01 ± 0.84 

-1.80 ± 0.36 

-1.01 ±0.66 

07III 

UVES 

HD 

167264 

10.4557 

-01.7408 

5.35 

0 

0.2 

1.60 ± 0.57 

-1.51 ± 0.41 

0.1 ± 0.45 

B0.5Ia/Iab 

UVES 

HD 

169454 

17.5385 

-00.6697 

6.65 

0.74 

0.93 

2.21 ± 0.87 

-4.49 ± 0.61 

-1.13 ± 0.72 

B1 la Ce 

UVES 

HD 

170235 

7.9523 

-06.7282 

6.59 

0.07 

0.28 

-1.00 ± 0.62 

-6.37 ± 0.39 

1.07 ± 0.61 

B2Vnne 

UVES 

HD 

188294 

32.6528 

-17.7667 

6.44 

-0.02 

0.09 

2.13 ± 0.49 

-30.77 ± 0.46 

6.42 ± 0.58 

B8V 

UVES 

HD 

199728 

28.537 

-36.5172 

6.25 

-0.13 

- 

10.07 ± 0.44 

-20.09 ± 0.36 

6.09 ± 0.42 

B9p 

UVES 

HD 

221507 

355.0398 

-70.3628 

4.37 

-0.08 

- 

95.97 ± 0.10 

38.29 ± 0.12 

18.74 ± 0.15 

B9.5IVmnp ... 

UVES 

HD 

223640 

60.5398 

-73.9377 

5.18 

-0.12 

- 

26.82 ± 0.38 

-4.27 ± 0.3 

10.23 ± 0.31 

Ap... 

UVES 

HD 

91316 

234.8870 

-1-52.7674 

3.870 

-0.15 

0.04 

-5.93 ±0.20 

-3.40 ±0.11 

0.60 ± 0.18 

Bllab 

McD 

HD 

114330 

311.4176 

-1-57.0274 

4.381 

0 

-0.01 

-36.28 ±0.93 

-31.22 ±0.69 

10.33 ± 1.09 

AlIVs 

McD 

HD 

141637 

346.0979 

-1-21.7059 

4.638 

-0.072 

0.15 

-14.20 ±0.33 

-25.12 ±0.26 

6.59 ± 0.27 

B1.5Vn 

McD 

HD 

143018 

347.2148 

-1-20.2305 

2.89 

-0.19 

0.04 

-11.42 ±0.78 

-26.83 ±0.74 

5.57 ± 0.64 

BIV 

McD 

HD 

143275 

350.0969 

-1-22.4904 

2.291 

-0.086 

0.14 

-10.21 ±1.01 

-35.41 ±0.71 

6.64 ± 0.89 

B0.2IVe 

McD 

HD 

144217 

353.1929 

-1-23.5996 

2.62 

-0.07 

0.17 

-5.20 ±0.92 

-24.04 ±0.64 

8.07 ± 0.78 

B0.5V 

McD 

HD 

144470 

352.7497 

-1-22.7730 

3.946 

-0.046 

0.18 

-8.98 ±0.23 

-23.48 ±0.16 

6.92 ± 0.26 

BIV 

McD 

HD 

147084 

352.3279 

-1-18.0503 

4.55 

0.84 

0.24 

-4.32 ±0.52 

-14.15 ±0.40 

3.71 ± 0.54 

A4H/HI 

McD 

HD 

147165 

351.3130 

-1-16.9989 

2.912 

0.097 

0.31 

-10.6 ±0.78 

-16.28 ±0.43 

4.68 ± 0.60 

B2IH 

McD 

HD 

148184 

357.9328 

-1-20.6766 

4.42 

0.28 

0.49 

-5.41 ±0.31 

-21.12 ±0.24 

6.21 ± 0.23 

B1.5Ve 

McD 

HD 

149438 

351.5348 

-1-12.8081 

2.814 

-0.201 

0.49 

-9.89 ±0.61 

-22.83 ±0.55 

6.88 ± 0.53 

BOV 

McD 

HD 

149757 

006.2812 

-1-23.5877 

2.578 

0.017 

0.30 

15.26 ±0.26 

24.79 ±0.22 

8.91 ± 0.20 

09.5Vnn 

McD 

HD 

152614 

028.7341 

-1-30.6605 

4.38 

-0.08 

0.03 

-53.8 ±0.22 

-34.04 ±0.19 

13.3 ± 0.22 

B8V 

McD 

HD 

159561 

035.8937 

-1-22.5681 

2.1 

0.13 

0.03 

108.07 ±0.87 

-221.57±0.80 

67.13 ± 1.06 

A5III 

McD 

HD 

184915 

031.7709 

-13.2866 

4.960 

-0.03 

0.19 

1.63 ±0.20 

-2.65 ±0.10 

1.94 ± 0.20 

B0.5HI 

McD 

HD 

186882 

078.7096 

-1-10.2430 

2.900 

-0.019 

0.02 

44.07 ±0.46 

48.66 ±0.49 

19.77 ± 0.48 

B9.5HI 

McD 

HD 

197345 

084.2847 

-1-01.9975 

1.25 

0.09 

0.02 

2.01 ±0.34 

1.85 ±0.27 

2.31 ± 0.32 

A2Iae 

McD 

HD 

198183 

078.0844 

-04.3378 

4.75 

-0.12 

0.03 

14.71 ±0.32 

-8.96 ±0.49 

4.24 ± 0.43 

B6IV 

McD 

HD 

200120 

088.0296 

-1-00.9707 

4.74 

-0.05 

0.17 

8.12 ±0.42 

2.15 ±0.50 

2.3 ± 0.42 

B1.5Vnne 

McD 

HD 

202904 

080.9767 

-10.0526 

4.43 

-0.11 

0.10 

10.03 ±0.33 

6.49 ±0.38 

5.08 ± 0.55 

B2Ve 

McD 

HD 

212571 

066.0067 

-44.7394 

4.794 

-0.168 

0.06 

17.83 ±0.27 

2.41 ±0.20 

4.17 ± 0.28 

BlVe 

McD 


in Table Most distances were calculated from the par¬ 
allax recorded in van Leeuwen (2007) as these tended to 
be the most reliable. The distance estimates from IHunterl 


et al. (2006) were calculated using the spectroscopic par¬ 


allax methodology given in Diplas & Savage (19941. The 


distances, along with the proper motions of the stars (van 


Leeuwen 2007) were used to calculate the transverse dis¬ 


placement on the sky during the interval between epochs 
by trigonometry, to allow us to set some limits on the scale 
over which any variation in the ISM occurred. We must note 


that there is no way to account for the motion of the clouds 
themselves which could act to increase or decrease this esti¬ 
mate on the scale of the structure in question, for example, 
if we assume the cloud itself to have a transverse velocity of 
Skms”'^, then in 6 yr it will have moved « 6 au, and de¬ 
pending on the direction of this motion, it could increase or 
decrease our transverse distance estimates calculated using 
the proper motion of the background star. We can therefore 
only say that the distance the star has moved across the sky 
is an indicator of the true scale of the ISM structure. Tabled 
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Table 1. Interstellar absorption-lines investigated in this paper, 
given in order of increasing wavelength. The wavelengths Aair 
and oscillator strengths / are taken f rom|Morton| l |2003| l for the 
atomic species, [Black Sz van Dishoecl^ < |1988[l for CH (mean of two 
unresolved transitions) and for CN, and|Gredel, van Dishoeck &| 
|Black| | |1993] for CH+. 


Trans. 

'^air 

(A) 

/ 

Nal 

3302.369 

0.009 

Nal 

3302.978 

0.005 

Till 

3383.759 

0.358 

Fel 

3719.935 

0.041 

Til 

3646.286 

0.236 

Can K 

3933.661 

0.627 

Can H 

3968.470 

0.312 

CN 

3873.999 

0.023 

CN 

3874.607 

0.034 

CH+ 

3957.692 

0.003 

Ki 

4044.142 

0.006 

Cal 

4226.728 

1.770 

CH+ 

4232.548 

0.006 

CH 

4300.313 

0.005 

Nal D 2 

5889.951 

0.641 

Nal Di 

5895.924 

0.320 


lists our final sample, whilst Fig.[2shows a histogram of the 
transverse distance moved for the 104 stars between the two 
epochs. 


2.5 Column Densities 


For sightlines where variations were detected, column den¬ 


sities were calculated using the vapid suite (Howarth et al. 
|2002[ ). This assumes Gaussian line of sight velocity distribu¬ 
tions of pure absorbers in each cloud ( Str6mgren| 1948 1, and 
so can obtain values for column density (N), central veloc¬ 
ity (v) and velocity dispersion (6) for each component of a 
specified absorption profile by fitting multiple Voigt profiles. 
Even at the high spectral resolution of the current study it 
is likely that most components will be unresolved at least in 
the UVES data set. WHK investigated the statistics of line 
widths and component separations at a significantly higher 
resolution than our UVES data. They found a median &(Na 
l) of 0.73 kms“^and median separation between adjacent 
components of about 2 kms“^(their figs. 5 and 7). Even at 
such high resolutions, these values are expected to be over 
estimates of the true values, and so the likelihood of unre¬ 
solved structure in our UVES data is very high. 

For Na I, the Di and D2 lines were fitted simultaneously, 
as were the H and K lines of Ca II. As very few velocity 
components showed variation between epochs, these stable 
components were measured in epoch one and fixed when 
fitting epoch 2, allowing us to accurately quantify the change 
in column densities in the varying components. 

Column densities for all varying sightlines were recal¬ 
culated, changing the instrumental resolution by 10 % and 
fixing the velocity values within vapid, allowing the 6- value 
and column density to vary freely. In most cases the differ¬ 
ences derived due to the change in resolution were negligible, 
and well within the errors presented for the column densi¬ 


ties. The only sightlines where larger variations due to this 
change in resolution were found, were in the particularly 
narrow components. However, even in these cases the varia¬ 
tion caused by changing the instrumental fo-value was small 
when compared with the variation seen between epochs. 


3 RESULTS 


In this section the results of the search for time variability 
in the twin-epoch spectra are discussed. We note that O- 
and B-type stars have high multiplicity (Mason et al.|2009 


Sana et al. 20131, and so great care must be taken not to 


confuse variations in stellar or circumstellar features with 
interstellar ones. Each pair of spectra and their residuals in 
the transitions listed in Table were searched by eye for ev¬ 
idence of differences between epochs. Only the strong lines 
of Nal D, Can H&K and the weaker Cal showed convinc¬ 
ing evidence for significant interstellar variation. The Nal 
D lines and Call H&K transitions each form a doublet, so 
putative variations were required to be present in both com¬ 
ponents of the doublet before they were accepted as real. For 
Cal, the variations are always seen at the same velocity as 
corresponding variations in Nal or Can which provides con¬ 
fidence in their reality. 

Of the total 83 twin-epoch UVES spectra analysed, 80% 
show no significant difference between the two epochs in any 
species. 2% show variation caused by eclipsing companions, 
12% show changes that are probably due to circumstellar 
material, and 6% show possible time variation in one or more 
velocity component caused by interstellar absorption - typ¬ 
ically there are > 5 velocity components per sightline. The 
sightlines where stellar or circumstellar variation are visible 
could also exhibit interstellar variations but they would be 
hidden due to the more prominent stellar-circumstellar pro¬ 
file changes. For the McDonald data, none of the sightlines 
show any substantial or convincing interstellar time vari¬ 
ation. The following subsections discuss each of the time- 
varying spectra. 


3.1 Sightlines that show time variation in stellar 
or circumstellar spectra 


Several of our sightlines show variations in the stellar spectra 
or circumstellar environment. Four such stars (HD 125823, 
HD 137509, HD 199728 and HD 22364) are classified as 


magnetic chemically peculiar stars ([Kochukhov & Bagnulo| 


20061, while HD 145792 is classified as a chemically pecu¬ 
liar, helium strong star ( Catanzaro|2008 1, and so variations 
in the spectra between epochs are to be expected. Some of 
the other stars have been grouped by their spectral classifi¬ 
cation below, explaining their spectral variations. 


5.1.1 Binary systems HD 123515 and HD 45116 

HD 123515 is a short-period double-lined spectroscopic bi¬ 
nary pulsating star l Aerts et al.||199^ . Figure]^ shows CH 
(43OOA) and Nal D2 (5895A) spectra towards the object at 
the two epochs. HD 47116 is also a star in a double system 
as can be seen from its spectra in Fig. 
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Table 3. Observation dates, derived distances d and transverse displacements for all stars in our sample. For some objects distances 
have been taken from the literature. References for these are given below the table. 


Star 

Distance (d) 

(pc) 

ref. 

Epoch 1 

date of observation 

epoch 2 

date of observation 

Interval 

(yr) 

Traversed distance 
(au) 

HD 10840 

180 ± 11.7 

P 

14-09-2002 

03-06-2008 

5.7 

28.3 ± 1.9 

HD 29138 

3530 ± 1059 

s 

21-01-2002 

20-06-2008 

6.4 

210 ± 63.9 

HD 30677 

1124 ± 669 

p 

28-02-2002 

26-07-2008 

6.4 

36.6 ± 22.0 

HD 33328 

249 ± 11 

p 

15-10-2002 

04-04-2008 

5.5 

2.70 ± 0.2 

HD 43285 

200 ± 17 

p 

23-02-2002 

28-08-2008 

6.5 

22.8 ± 2.0 

HD 45725 

184 ± 55 

s 

13-03-2001 

04-04-2008 

7.1 

11.2 ± 3.6 

HD 46185 

380 ± 80 

p 

16-02-2002 

14-08-2008 

6.5 

10.4 ± 2.6 

HD 46462 

341 ± 61 

p 

25-02-2002 

14-08-2008 

6.5 

28.5 ± 5.2 

HD 47116 

654 ± 261 

p 

25-02-2002 

02-09-2008 

6.5 

23.2 ± 9.5 

HD 48099 

855 ± 300 

p 

09-04-2001 

19-09-2008 

7.5 

16.8 ± 6.2 

HD 49131 

581 ± 183 

p 

02-10-2002 

18-08-2008 

5.9 

15.5 ± 5.2 

HD 49333 

242 ± 30 

p 

13-09-2002 

18-08-2008 

5.9 

22.1 ± 2.8 

HD 50896 

1389 ± 1196 

p 

20-01-2002 

17-09-2008 

6.7 

53.2 ± 46.0 

HD 51876 

333 ± 66 

p 

14-02-2003 

13-09-2008 

5.6 

16.2 ± 3.3 

HD 52918 

373 ± 31 

p 

26-10-2002 

04-04-2008 

5.4 

11.2 ± 1.1 

HD 58343 

287 ± 48 

p 

13-03-2001 

04-04-2008 

7.1 

42.2 ± 7.1 

HD 58377 

448 ± 134 

s 

14-02-2003 

03-09-2008 

5.6 

22.7 ± 6.9 

HD 58978 

383 ± 48 

p 

03-10-2002 

21-09-2008 

6.0 

20.9 ± 2.7 

HD 60498 

239 ± 19 

p 

15-02-2003 

19-09-2008 

5.6 

12.8 ± 1.1 

HD 61429 

190 ± 9.4 

p 

29-11-2001 

04-04-2008 

6.4 

9.66 ± 0.5 

HD 62714 

244 ± 24 

p 

25-02-2002 

17-09-2008 

6.6 

25.6 ± 2.7 

HD 64972 

2858 ± 857 

s 

27-11-2002 

17-04-2008 

5.4 

175 ± 56.9 

HD 67536 

368 ± 30 

p 

26-10-2002 

17-04-2008 

5.5 

25.4 ± 2.1 

HD 68761 

613 ± 128 

p 

25-02-2002 

17-04-2008 

6.1 

27.3 ± 5.8 

HD 72067 

439 ± 94 

p 

17-11-2001 

17-04-2008 

6.4 

27.8 ± 6.2 

HD 74966 

446 ± 104 

p 

11-12-2002 

17-04-2008 

5.4 

29.2 ± 6.8 

HD 76131 

420 ± 51 

p 

15-02-2003 

17-04-2008 

5.2 

38.3 ± 4.7 

HD 76341 

2365 ± 709 

s 

24-02-2002 

17-04-2008 

6.1 

90.6 ± 27.8 

HD 83312 

518 ± 107 

p 

15-01-2003 

17-04-2008 

5.3 

41.7 ± 8.7 

HD 88661 

391 ± 56 

p 

08-02-2001 

17-04-2008 

7.2 

35.3 ± 5.2 

HD 89587 

900 ± 270 

s 

30-11-2001 

17-04-2008 

6.4 

73.2 ± 22.1 

HD 90177 

255 ± 76 

s 

08-10-2002 

07-06-2008 

5.7 

10.5 ± 3.4 

HD 90264 

123 ± 2.7 

p 

26-02-2002 

03-04-2008 

6.1 

18.9 ± 0.4 

HD 90882 

99 ± 2.2 

p 

24-11-2001 

03-04-2008 

6.4 

31.8 ± 0.7 

HD 92740 

2600 ± - 

b 

02-01-2003 

17-04-2008 

5.3 

113 ± - 

HD 94910 

6000 ± 1000 

c 

11-01-2003 

08-06-2008 

5.4 

207 ± 43.6 

HD 94963 

3257 ± 977 

s 

25-02-2002 

17-04-2008 

6.1 

168 ± 56.7 

HD 100841 

129 ± 5.6 

p 

23-06-2001 

17-04-2008 

6.8 

30.0 ± 1.3 

HD 105056 

5184 ± - 

b 

23-02-2002 

26-06-2008 

6.3 

149 ± - 

HD 105071 

3482 ± 1044 

s 

26-02-2002 

17-04-2008 

6.1 

124 ± 38.4 

HD 106068 

2824 ± 847 

s 

20-07-2001 

17-04-2008 

6.7 

97.1 ± 29.7 

HD 109867 

3264 ± 979 

s 

26-02-2002 

17-04-2008 

6.1 

93.6 ± 28.7 

HD 112272 

714 ± 270 

p 

26-02-2002 

17-04-2008 

6.1 

53.1 ± 20.2 

HD 112842 

2126 ± 638 

s 

15-02-2003 

17-04-2008 

5.2 

47.8 ± 15.2 

HD 113904 

2817 ± 845 

s 

22-07-2001 

17-04-2008 

6.7 

90.9 ± 28.4 

HD 115363 

3282 ± 984 

s 

26-02-2002 

26-06-2008 

6.3 

177 ± 54.1 

HD 115842 

2157 ± 647 

s 

25-02-2002 

17-04-2008 

6.1 

65.9 ± 20.3 

HD 120709 

105 ± 10 

p 

24-06-2001 

03-04-2008 

6.8 

31.9 ± 3.0 

HD 123515 

251 ± 25 

p 

05-08-2002 

17-04-2008 

5.7 

49.1 ± 4.8 

HD 125823 

140 ± 3.1 

p 

13-03-2001 

03-04-2008 

7.1 

32.3 ± 0.7 

HD 133518 

452 ± 90 

p 

13-03-2001 

17-04-2008 

7.1 

41.2 ± 8.4 

HD 136239 

2350 ± 705 

s 

22-02-2002 

17-04-2008 

6.2 

88.4 ± 30.0 

HD 137509 

196 ± 15 

p 

11-07-2001 

17-04-2008 

6.8 

28.8 ± 2.2 

HD 137753 

303 ± 39 

p 

14-02-2003 

17-04-2008 

5.2 

18.6 ± 2.4 

HD 142301 

158 ± 11 

p 

13-03-2001 

09-04-2008 

7.1 

32.1 ± 2.3 

HD 142758 

1316 ± 883 

p 

26-02-2002 

17-04-2008 

6.1 

61.2 ± 41.3 


Distance reference codes; p - distances calculated from the parallax recorded in 

van Leeuwen| ||2007|). 

s - distances calculated spectroscopically ||Hunter 3t al.|2006 

1. 


b - |Diplas fc Sava"^ | |199^ , c - |Hoekzema, Lamers Sz van Genderen| | |1993[ |, d - |Conti &: Vac^ | |1990| l 
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Table 3. continued 


Star 

Distance (d) 

(pc) 

ref. 

Epoch 1 

date of observation 

epoch 2 

date of observation 

Interval 

(yr) 

Traversed distance 
(au) 

HD 142983 

143 ± 4.9 

P 

11-07-2001 

03-04-2008 

6.7 

20.1 ± 0.7 

HD 143448 

649 ± 245 

P 

26-02-2002 

17-04-2008 

6.1 

18.9 ± 7.4 

HD 145482 

147 ± 3.5 

P 

13-08-2002 

09-04-2008 

5.7 

21.7 ± 0.5 

HD 145792 

144 ± 13 

P 

06-10-2002 

09-04-2008 

5.5 

17.5 ± 1.7 

HD 145842 

118 ± 3.6 

P 

26-02-2002 

03-04-2008 

6.1 

40.8 ± 1.3 

HD 148184 

161 ± 6.0 

P 

26-02-2002 

03-04-2008 

6.1 

21.4 ± 0.8 

HD 148379 

559 ± 218 

P 

26-02-2002 

09-04-2008 

6.1 

16.4 ± 7.4 

HD 148688 

833 ± 229 

P 

14-09-2002 

09-04-2008 

5.6 

21.2 ± 6.1 

HD 148937 

426 ± 143 

P 

26-02-2002 

09-04-2008 

6.1 

10.1 ± 4.0 

HD 152003 

2687 ± 806 

S 

26-02-2002 

09-04-2008 

6.1 

58.1 ± 24.0 

HD 152235 

952 ± 517 

P 

03-10-2002 

09-04-2008 

5.5 

18.2 ± 10.1 

HD 154811 

427 ± 126 

P 

03-09-2002 

09-04-2008 

5.6 

9.63 ± 3.0 

HD 154873 

463 ± 154 

P 

07-10-2002 

09-04-2008 

5.5 

12.9 ± 4.5 

HD 155416 

1278 ± 383 

s 

27-07-2001 

09-04-2008 

6.7 

20.0 ± 6.4 

HD 155806 

1064 ± 319 

s 

11-07-2001 

09-04-2008 

6.8 

22.9 ± 7.0 

HD 156385 

1300 ± - 

d 

14-10-2002 

17-04-2008 

5.5 

30.6 ± - 

HD 157038 

1444 ± 433 

s 

22-09-2002 

03-06-2008 

5.7 

11.0 ± 5.4 

HD 163745 

380 ± 58 

p 

19-06-2001 

27-05-2008 

6.9 

10.6 ± 1.8 

HD 163758 

4103 ± 1231 

s 

27-07-2001 

09-04-2008 

6.7 

160 ± 58.9 

HD 163800 

1621 ± 486 

s 

29-09-2002 

23-06-2008 

5.7 

16.7 ± 6.0 

HD 167264 

2109 ± 632 

s 

12-07-2001 

27-05-2008 

6.9 

31.9 ± 12.0 

HD 169454 

1359 ± 407 

s 

13-09-2002 

23-06-2008 

5.8 

39.3 ± 12.9 

HD 170235 

501 ± 150 

s 

10-07-2001 

19-04-2008 

6.8 

21.9 ± 7.0 

HD 188294 

156 ± 14 

p 

12-07-2001 

23-04-2008 

6.8 

32.6 ± 3.0 

HD 199728 

164 ± 11 

p 

29-11-2002 

23-06-2008 

5.6 

20.6 ± 1.5 

HD 221507 

53.4 ± 0.4 

p 

09-06-2002 

27-05-2008 

6.0 

32.9 ± 0.3 

HD 223640 

97.8 ± 3 

p 

13-08-2002 

03-06-2008 

5.8 

15.4 ± 0.5 

HD 091316 

1346 ± 403 

s 

- 

29-05-2008 

20.0 * 

184 ± 55.3 

HD 114330 

97 ± 10 

p 

27-05-2008 

04-07-2014 

6.1 

28.3 ± 3.0 

HD 141637 

152 ± 6.2 

p 

23-06-1988 

28-05-2008 

19.9 

167 ± 50.1 

HD 143018 

180 ± 21 

p 

21-06-1988 

28-05-2008 

19.9 

104 ± 30.9 

HD 143275 

151 ± 20 

p 

31-05-2008 

06-07-2014 

6.1 

38.9 ± 11.7 

HD 144217 

124 ± 12 

p 

- 

27-05-2008 

20.0 * 

77.2 ± 23.2 

HD 144470 

145 ± 5.4 

p 

01-05-2008 

06-07-2014 

6.2 

37.2 ± 11.0 

HD 147084 

270 ± 39 

p 

28-05-2008 

04-07-2014 

6.1 

24.3 ± 3.6 

HD 147165 

214 ± 27 

p 

28-05-2008 

04-07-2014 

6.1 

22.1 ± 6.6 

HD 148184 

161 ± 6.0 

p 

30-05-2008 

04-07-2014 

6.1 

21.4 ± 0.8 

HD 149438 

145 ± 11 

p 

21-06-1988 

30-05-2008 

20.0 

72.0 ± 5.7 

HD 149757 

112 ± 2.5 

p 

10-06-1905 

28-05-2008 

20.0 * 

116 ± 34.4 

HD 152614 

75 ± 1.2 

p 

30-05-2008 

06-07-2014 

6.1 

29.2 ± 0.5 

HD 159561 

15 ± 0.2 

p 

28-05-2008 

06-07-2014 

6.1 

22.6 ± 22.6 

HD 184915 

515 ± 53 

p 

27-05-2008 

04-07-2014 

6.1 

12.4 ± 3.7 

HD 186882 

51 ± 1.2 

p 

- 

27-05-2008 

20.0 * 

49.9 ± 14.4 

HD 197345 

433 ± 60 

p 

29-05-2008 

06-07-2014 

6.1 

7.22 ± 1.29 

HD 198183 

236 ± 24 

p 

29-05-2008 

06-07-2014 

6.1 

24.8 ± 2.6 

HD 200120 

435 ± 79 

p 

23-06-1988 

31-05-2008 

20.0 

49.6 ± 15.1 

HD 202904 

197 ± 21 

p 

11-10-1989 

30-05-2008 

18.6 

53.2 ± 15.9 

HD 212571 

240 ± 16 

p 

31-05-2008 

06-07-2014 

6.1 

45.8 ± 13.7 


* epoch interval presumed to be 20 yr as no specific observation date available for epoch 1 (observed between 1987 and 1989) 


Distance reference codes: p - distances calculated from the parallax recorded in 

van Leeuwen| ||2007|). 

s - distances calculated spectroscopically ifHunter 3t al.|2006 

. 


b -|Diplas fc Savag^ ||1994ll, c -|Hoekzema, Lamers &: van Genderen] ||1993[l, d - [Conti fe Vac^(|1990[l 


3.1.2 B(e) stars HD 58978, HD 83312, HD 90177, 

HD 94910 and HD I 48 I 84 

These five objects, whose spectra show significant variation, 
are in fact B(e) stars. Two examples are HD 58978 and 
HD 83312 whose twin-epoch spectra are shown in Figs. 
and respectively. HD 58978 (FY CMa) is a well-studied 
object and was found to be a B(e)-|-sdO binary with a period 


of 37.26±0.03 days by Rivinius et al. (20041. It shows strong 
shell phases and ionization in the outer parts of the circum- 
stellar disc. Less studied is the emission-line star HD 83312 
|Henize||l976 1 which is part of the catalogue of young run¬ 


away Hipparcos stars within 3 kpc of the Sun (Tetzlaff, 
Neuhauser &: Hohle||20lf I. The spectra of a typical Be star 


often shows clumpy winds with discrete accelerating compo- 
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Figure 2. Overlay of epoch 1 (solid red line) and epoch 2 (dashed 
blue line) spectra for HD 123515, showing the areas of spectrum 
in the range of the CH (above) and Na D 2 lines (below). The 
line profiles and varations visible are stellar lines varying due to 
stellar binarity, making interstellar variation almost impossible to 
detect. Many telluric lines are also seen in the Na D 2 spectra. 
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Figure 3. Overlay of epoch 1 (solid red line) and epoch 2 (dashed 
blue line) spectra for HD 47116, showing variations in line profile 
of the CH (above) and Na D 2 lines (below), due to stellar binarity. 


HD-58g78 



Figure 4. Overlay of epoch 1 (solid red line) and epoch 2 (dashed 
blue line) spectra for HD 58978, showing variations in line profile 
of the Call K (above) and Can H lines (below) due to B(e) star 
variability. 
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Figure 5. Overlay of epoch 1 (solid red line) and epoch 2 (dashed 
blue line) spectra for HD 83312, showing variations in line profile 
of the Ca II K (above) and Ca II H lines (below) due to B(e) star 
variability. 


nents ( Stefl et al.||2012 and references therein). HD 94910, 
HD 90177 and HD 148184 have also been widely studied, 
(e.g. Burbidge 1952[[Henize|1957 1 and show the similar spec¬ 
tral characteristics of B(e) stars. Any ISM lines have profiles 
so compromised by the stellar variation in such sightlines, 
that we therefore classify the changes in these spectra as cir- 
cumstellar in nature, and their ISM lines are not discussed 
further. 


3.2 Possible time-variable spectra in the 
interstellar medium 

In total fix targets (6%) from the UVES sample show appar¬ 
ent temporal variations in interstellar lines, which are un¬ 
likely to be explained by stellar changes. A summary of all 
the sightlines where possible variation is seen can be found 
in Table 1^ and are discussed individually below. Although it 
is possible that some variations may be caused by slight res¬ 
olution differences between epochs (see Sect. [T5| ), had this 
been the case it would be expected that variation would be 
much more frequently observed throughout the data set. As 
80 % show no variation we are convinced that these five 
sightlines show genuine evidence of fine-scale structure. 


3.2.1 HD 112842 - possible Nal D and Call H & K 
variations 

Possible variations can be seen in intermediate velocity gas 
in the profiles of both the Nal D and in Call H & K 
lines in this sight line (Fig. and Fig. [^, with the com¬ 
ponents at —40kms“^ stronger in epoch 1 in all cases. The 
proper motion of this star is 9.2 an yr“^, and there are 5.2 
yr between epochs 1 and 2. Hence the star has travelled 
~48 an perpendicular to the line of sight between epochs. 
The column density of the —40kms“^ components decrease 
from log Af(Na I (cm“^))=11.45 ± 0.01 to 11.36 ± 0.02, 
logiV(Call (cm“^))=11.72 ± 0.12 to 11.33 ± 0.11, or a 
change of 19 percent in Na I and a factor of 2.5 in Ca ll. 
No variations are seen in other species. For this sightline we 
have obtained 3rd epoch observations on 29 November 2013, 
which are also plotted in Fig. and Fig.[^ It can clearly be 
seen in Call H & K that there is no variation at all in the 
5.6 yr and 51.9 an between epoch 2 and 3, and so residuals 
plotted are between epoch 1 and 2. For Nal D, resolution 
changes are visible between 2008 and 2013, making it im¬ 
possible to quantify if there is any real variation in the 3rd 
epoch, and so the residuals plotted are once again, only be¬ 
tween the earlier two epochs. To verify that the variation 
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is not solely due to differences in the resolution between 
epochs, Gaussian smoothing was applied to the data in an 
effort to make the effective instrnmental resolution the same 
in all epochs. It was fonnd that changes in resolntion could 
not force all components to match in all epochs, implying 
that the variation is real. Ca l may also show a decrease 
between epoch 1 and 2, althongh the quality of the epoch 
2 data is not optimal at this wavelength to see convincing 
variation and the Ca l lines are weak. 

HD 112842 is at a distance of ~2.1±0.5 kpc and lies 
in the Centaurus OBI stellar association. Observations of 
Hi towards this region have revealed a large expanding 
shell (GSH 305+04-26) with a velocity range of —33 to 
— 17kms“^ and estimated distance of 2.5+0.9 kpc and di¬ 
mensions of 440 X 270 pc ( Corti, Arnal fc Orellan^|2012 1. 
As HD 112842 is predicted to be closer than the Hi shell, 
this implies there is fine-scale structnre in the normal dif¬ 
fuse ISM. However, the errors on the distances are large and 
so there remains the possibility that the variation seen is 
cansed by part of the shell. Nal absorption is detected in 
the velocity range, although the variable component is more 
blue-shifted than expected for the shell. 


3.2.2 HD 113904 - Possible Nal D variations 

HD 113904 {6 Mus) is a WR star with a binary compan¬ 
ion (WR 48 WC5+06/7V) and possibly has an additional 
O9.5/B0Iab-type supergiant companion ( [Sugawara, TsuboI| 
&: Maeda||2008 1 , located in the same area of the sky as the 
Centaurus OBI association. However the proper motion of 
the star is in the opposite direction to the members of this 


association. Hence, Corti, Arnal & Orellana (20121 deem it is 
unlikely to be part of Cen OBl. |Stupar, Parker fc Filipovic| 
(20101 give the distance as between 0.9 and 1.3 kpc, and 


from Hq imaging postulate that the nebulosity seen near 
to the star is not a WR shell but in fact a possible optical 
supernova remnant and complex of Hll regions. 

The bluest-shifted component, (see Fig. |^, at 
—60kms“^, shows possible variations in intermediate ve¬ 
locity gas in both Nal lines in this sightline, with stronger 
absorption at epoch 2. As the proper motion of this star is 
13.48 auyr“^, and there are 6.7 yr between observations, 
the transverse displacement is ~91 au. Although these vari¬ 
ations are small - logA'(Nal) for this component increases 
from 11.52 ± 0.01 dex to 11.56 ± 0.01 dex (10 percent) - 
they can be seen in the residual plots. For this sightline, we 
have obtained 3rd epoch observations from 29 November, 
2013 which are also plotted in Fig. No variation is visible 
over the 75.8 au traversed between 2008 and 2013, so the 
residuals plotted are only between the earlier two epochs. 
There are some instrumental effects acting upon these spec¬ 
tra, as we can see the resolution is not identical for all three 
epochs. However, the variation we do see is strong enough 
to be identified despite these slight differences in resolution. 


3.2.3 HD 163758 - Possible Nal D variations 

HD 163758 is a WR star of spectral type 06.5Iaf at a dis¬ 
tance of 4.1 ± 1.2 kpc. For this sightline we were able to 
acquire epoch 3 data, observed on 19 September 2013. In 
Fig.i we can clearly see variations in two of the compo¬ 
nents of both Nai D lines, at —26kms“^ and -1-19 kms“^ 


HD 112842 



Figure 6. Overlay of epoch 1, 2002 (solid red line) and epoch 
2, 2008 ( dash-dotted blue line) and epoch 3, 2013 (dotted black 
line) spectra for HD 112842, showing variations at —40kms“^in 
the profiles of the Can K (top plots) and Call H (bottome plots) 
lines. The residuals between epochs 1 and 2 are plotted below 
each line. 


between all three epochs. Both components strengthen with 
time, that at +19 km s“^ showing the largest increase from 

2002 to 2008 (logiV(Na (cm"H)= 12-08 ± 2.0 to 12.63 

± 2.0, a change of a factor of 3.55) with very little vari¬ 
ation between 2008 and 2013. The particularly small b- 
values (of « 0.5 kms“^) of these components, mean they 
are unresolved and therefore uncertain, leading to large 
errors in the column densities, which should be regarded 
as approximate estimates. The component at —26kms“^ 
strengthens by a small amount between 2002 and 2008, 
logN(Na (cm"H)=ll-87 ± 0.02 to 11.90 ± 0.01, a 7% 
change, and then by a much larger degree between 2008 and 
2013, to 11.99 ± 0.01, a 32% change. These variations can 
very clearly be seen in the residual plots of each Nal D line. 
Between epoch 1 and 2 the star moves 160 au in 6.7 yr and 
in the 5.4 yr between epoch 2 and 3 travels 130 au. We note, 
there are some slight differences in resolution between the 
three epochs due to instrumental effects, as can be seen, for 
example at 0 km s“^, possibly bringing the weaker variations 
into question. As the strongly varying component is very 
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HD 112842 



Figure 7. Overlay of epoch 1, 2002 (solid red line) and epoch 2, 
2008 ( dash-dotted blue line) and epoch 3, 2013 (dotted black line) 
spectra for HD 112842, showing variations at —40kms“^in the 
profiles of the Na D lines (D 2 above and Di below). The residuals 
between epochs 1 and 2 are plotted below each line. Telluric lines 
have been removed. 


HD 113904 



Figure 8. Overlay of epoch 1, 2002 (solid red line) and epoch 2, 
2008 ( dash-dotted blue line) and epoch 3, 2013 (dotted black line) 
spectra for HD 113904, showing variations at —60kms“^in the 
profiles of the Na D lines (D 2 above and Di below). The residuals 
between epochs 1 and 2 are plotted below each line. Telluric lines 
have been removed. 


narrow, the changes in resolution have more effect on this 
component. However, even varying the spectral resolutions 
by 10 percent cannot account for the stronger variations that 
we see, (see Sect. 2.51. 


3 . 2.4 tlD 50896 - possible Nai D variations 


In the sightline of the well-studied WR star HD 50896 
(WR 6) (see Fig. lOl, the component at 21.3kms“^ has 
strengthened in both Nai D lines by a factor of 2.3, from 
log Ar(Na (cm“D) = 12-50 ± 2.0 to 12.87 ± 2.0 in epoch 2, 
showing possible variation in the 6.7 yr between the two 
observations. These components have particularly small b- 
values (of « 0.5 kms“^), meaning they are unresolved and 
therefore uncertain, leading to large errors in the column 
densities, which should be regarded as approximate esti- 


the transverse displacement is only 53.2 au. These variations 
are small but tweaking the continuum fitting and resolution 


would only prove to increase the variations, and so we are 
sufficiently convinced that these variations are real. 


3.3 HD 29138 - Cal and Call variations 

The variations seen in the sightline of HD 29138 are the 
most dramatic in the whole sample, and can very definitely 
not be due to differences in spectral resolution, as changing 
the instrumental resolution of this sightline by 10 percent 
resulted in no changes in the column densities at all. 

HD 29138 is a Bllab star located far from the plane of 
the Milky Way {b = —30.5° and d — 3530pc). Significant 
variations in Ca l and both Ca ll H&K lines can be seen. All 
calcium lines show changes at —18.6kms“^ where a feature 
seems to diminish in all three lines in epoch 2 (Fig. |ll[ ). The 
change in the Call column density of this profile is 0.13 dex 
(or 35 percent), from 11.53 ± 0.05 to 11.40 ± 0.01 dex. As 
the Cal component disappears completely in epoch 2, we 
were only able to measure a column density for epoch 1. 
However, using the S/N of the spectra at both epochs (see 
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HD 163758 



Figure 9. Overlay of epoch 1, 2002 (solid red line) and epoch 
2, 2008 ( dash-dotted blue line) and epoch 3, 2013 (dotted black 
line) spectra for HD 163758, showing variations at —26kms“^ 
and 19kms“^ in the profiles of the Na D lines (D 2 above and 
Di below). The residuals between epochs are plotted below each 
line. 


HD 50896 



Figure 10. Overlay of epoch 1, 2002 (solid red line) and epoch 
2, 2008 ( dash-dotted blue line) and epoch 3, 2013 (dotted black 
line) spectra for HD 50896, showing variations at 21kms“^ in 
the profiles of the Nal D lines (D 2 above and Di below). The 
residuals between epochs are plotted below each line. 


Section 3.41 it was possible to calculate an upper limit for 


the column density of Ca l in epoch 2. The difference be¬ 
tween this upper limit and epoch 1 is 0.47 dex, i.e. the col¬ 
umn density decreases from 9.60 ± 0.02 to an upper limit of 
9.1 dex (change of a factor of 3.2). As the proper motion of 
this star is 32.8 an yr”'^ these changes in the profiles occur 
over 210 au, due to the 6.4 yr between epochs. Although the 
changes are relatively small, the fact that we see variations 
in three lines and two species at the same velocity indicates 
that they are real. We do not see obvious variation in Nal or 


Till (which is usually closely correlated with Call; (Welsh 
et al. 1997| Hunter et al. 20061) which could indicate cir- 


cumstellar variation (c.f. Kiefer et al. 20141. However, the 
Na D doublet is stronger than the Ca lines and so any small 
variation would be very difficult to see. epoch 3 data were 
obtained for this sightline on 19 September, 2013, 5.3 yr af¬ 
ter epoch 2 accounting for a transversed distance of 172 au 
but no further variation is observed. 

The epoch 1 data also shows an unresolved line corre¬ 
sponding to Fe I at 3859.9114 A, with central velocity of 


— 19.2 kms“^ and column density in the optically thin ap¬ 
proximation of logN (Fe l)= 11.66. The line disappears in 
the second epoch, with a 3cr upper limit to its column den¬ 
sity of logAf (Fe l)= 11.32. The epoch 3 data are too noisy 
to set firm limits on the Fe l. The Fe l line observed in epoch 
1 is close in velocity to that observed in Ca l (4226 A) which 
has V = -18.6 kms“^. We would expect to see a stronger Fe 
I feature at 3719.93 A due to its larger oscillator strength, 
however, no absorption is detected at this wavelength, cast¬ 
ing doubt on the reliability of the line observed at 3859.91 
A. Using the Nal Di line we have calculated an approxi¬ 
mate column density for Na l of logA^(Na (cm“^)) ~ 10.9 
dex and a Ca l abundance of logA^(Ca (cm~^)) = 9.6 ±0.1 
dex, yielding a relatively high Ca l/Na l ratio. If the observed 
change in line strength of Fe l is real, then the disappearance 
of both Ca I and Fe l, combined with the high Ca l/Na l ra¬ 
tio, suggest that the variable component is an example of the 


CaFe” clouds noted by Bondar et al. 120071 and discussed 


by Gnacinski & Krogulec 120081 and Welty, Simon & Hobbs 
| |2008p . The enhanced Ca i and Fe I in these cases could be 
due to some combination of dielectronic recombination and 
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milder than usual depletions at temperatures above « 5000 

K. 

3.4 Sightlines showing no time-variations 

The vast majority of the sightlines in our survey show little 
or no variation in line profile between the two epochs for 
the interstellar species isolated and examined. However at 
our spectral resolution (J?=80,000) there is always the pos¬ 
sibility that there are blends in the absorption-lines of our 
observations. 

We estimated upper limits to the changes, 
A log NEpoch 2 , in the column density of each interstel¬ 
lar component by using the the expression: 


Al0giVepoch2 = AlogAepochl^l^^^^, (1) 

10/rVepoch2 j 

where A log Aepochi is the error in the log of the column 
density of each component for epoch 1 from Papers I and 
V, S/Nepochi and S/Nepoch 2 are the signal to noise ratios 
from epoch 1 and epoch 2, respectively. The column density 
change limit has only been derived for Call K. 

Typical errors in the column density for epoch 1 are 
from 0.02 to 0.06 dex, being on average twice as large for 
epoch 2. The resulting Scr upper limits to changes in Call 
are thus typically 0.1 to 0.3 dex. 

None of our McDonald sightlines show very clear varia¬ 
tion. Two of this sample, namely HD 91316 and HD 149757, 
show some possible variation. However the difference in res¬ 
olution between epochs cast doubt on the reliability of such 
variations. As this variation is not as convincing as that of 
the UVES sightlines, they are not included in our final re¬ 
sults here, but some discussion of these sightlines is provided 
in Appendix [X] 


4 DISCUSSION 


In this section we discuss the main results of this survey 
and relate these to the idea of biased neutral formation in 


the ISM (Lauroesch & Meyer 20031. Excluding the sight¬ 


lines where we are likely seeing changes in absorption pro- 
Hles due to circumstellar material or companions, some 1% 
of the sightlines in Cal, 2% of those in Call, K and H, and 
4% in Nal D show evidence for time variation. None was de¬ 
tected in Till, Eel, CN, CH^ or Kl (4044A line only), most 
probably because the absorption features are too weak to see 
changes, due to depletion or smaller /-values, or in the case 
of Ti II due to a broader, smoother spatial distribution. The 
smallest scale structure over which a change in the ISM was 
tentatively observed was ~50 au, even though some 64 % of 
the stars in our sample moved a smaller distance between 
the two epochs. 

Optical and Hi observations of interstellar clouds im¬ 
ply gas densities several orders of magnitude higher than the 
background ISM (e.g. Diamond et al. 1989 Crawford|2002 


[Smith et al.|[^13[ ), so they are not in pressure equilibrium 
with the surrounding gas. Such clouds would quickly evap¬ 
orate, so must be continuously generated or confined by ad¬ 
ditional sources of pressure. [Heilesj ( |1997| ) invoked the idea 
of non-spherical geometry to solve this pressure problem. 


HD 29138 



Velocity (km/s) 

Figure 11. Overlay of epoch 1, 2002 (solid red line) and epoch 
2, 2008 ( dash-dotted blue line) and epoch 3, 2013 (dotted black 
line) spectra for HD 29138, showing variations at —19kms“^ in 
the profiles of calcium lines (Can H above, Can K middle, and 
Cal below). The residuals between epochs are plotted below each 
line. 
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Table 4. Table of species where we see possible variation in Na i D, Ca ii H and K and Ca i 


Sight line 

Species 

Component 

velocity 

Column density 
(1st epoch) 

Column density 
(2nd epoch) 

Column density 
variation (%) 

time 

(y) 

traversed 
distance (au) 

HD 112842 

Na I D 

-40.8 ± 0.20 

11.45 ± 0.01 

11.36 ± 0.02 

19 

5.2 

47.8 ± 15.2 


Ca II K, H 

-39.5 ± 0.20 

11.72 ± 0.12 

11.33 ± 0.11 

250 

5.2 

47.8 ± 15.2 

HD 113904 

Na I D 

-60.5 ± 0.09 

11.52 ± 0.01 

11.56 ± 0.01 

10 

6.7 

90.9 ± 28.4 

HD 163758 

Na I D 

-25.9 ± 0.05 

11.87 ± 0.02 

11.90 ± 0.01 

7 

6.7 

160.3 ± 58.9 


Na I Dt 

-26.0 ± 0.07 

11.90 ± 0.01 

11.99 ± 0.01 

32 

5.4 

130.2 ± 47.8 


Na I D 

18.3 ± 1.00 

12.08 ± 2.00 * 

12.63 ± 2.00* 

355 

6.7 

160.3 ± 58.9 

HD 50896 

Na I D 

21.5 ± 0.08 

12.50 ± 2.00* 

12.87 ± 2.00* 

233 

6.7 

53.2 ± 46.0 

HD 29138 

Ca I 

-18.6 ± 0.15 

9.60 ± 0.02 

9.1 * 

320 

6.4 

210.2 ± 63.9 


Ca II K, H 

-18.7 ± 0.11 

11.53 ± 0.05 

11.40 ± 0.01 

35 

6.4 

210.2 ± 63.9 


^ Variation seen between epoch 2 and 3. 

* Upper limit to column density value. 

★ Sightlines with particularly narrow 6-values, and therefore large errors 


Table 5. Table of velocities, 6-values and column densities of components towards stars where it was possible to calculate rig and lower 
limits for nu-, which are also listed in the table, nj/ estimates are uncertain, see Sect. [4T] for discussion of caveats. For sightlines where 
no variation was discovered, component velocities, 6-values and column densities were taken from [Hunter et al.| l |200'^ and Smoker & 
Ledoux (in preparation) 


Star 

Ca I 

V (kms“^) 

b (km 

s-l) 

logAf ( 

cm ^) 

Ca II 

V (kms~^) 

b (km s ^) 

logW (cm 

rie 
(cm 

-3 

) 

riH 

(cm“®) 


HD 163758 

- 26.7 ± 0.4 * 

3.4 

± 

0.7 

9.60 

± 

0.10 

-26.2 ±0.1 

4.0 

±0.1 

11.99 

±0.01 

0.27 ± 0.03 • 


1917 

± 

195 • 


-4.3 ± 0.1 

2.0 

± 

0.2 

10.13 

± 

0.02 

-3.3 ±0.1 

4.3 

±0.2 

12.48 

±0.02 

0.29 ± 0.01 


2086 

± 

0.5 


4.8 ± 0.2 

4.6 

± 

0.4 

10.13 

± 

0.02 

5.3 ±0.1 

4.8 

±0.2 

12.59 

±0.02 

0.23 ± 0.01 


1647 

± 

13 


24.3 ± 0.2 

2.9 

± 

0.4 

9.78 

± 

0.03 

25.4 ±0.1 

4.3 

±0.1 

11.91 

±0.01 

0.49 ± 0.03 


3497 

± 

186 

HD 29138 (1) 

-18.7 ± 0.1 * 

1.3 

± 

0.3 

9.60 

± 

0.02 

-18.6 

± 

0.1 

1.9 

± 

0.3 

11.50 

± 

0.05 

0.70 ± 0.10 * 


4515 

± 

375 • 

HD 29138 (2) 

-18.7 ± 0.1 * 

1.3 

± 

0.3 



9.10 

-18.5 

± 

0.1 

1.9 

t 


11.40 

± 

0.01 

^ 0.3 * 


- 




HD 30677 

12.1 

± 

0.2 

0.9 

± 

0.0 

9.93 

± 

0.08 

11.8 

± 

0.7 

4.0 

± 

0.3 

11.70 

± 

0.02 

1.12 

± 

0.17 


8006 

± 

1186 

HD 33328 

-11.5 

± 

0.1 

2.0 

± 

0.4 

9.71 

± 

0.08 

-11.8 

± 

0.7 

4.0 

± 

0.1 

11.86 

± 

0.01 

0.47 

± 

0.08 


3337 

± 

584 


0.8 

± 

0.5 

2.7 

± 

2.7 

9.25 

± 

0.30 

0.0 

± 

0.4 

2.2 

± 

0.1 

11.38 

± 

0.01 

0.49 

± 

0.46 


3495 

± 

3319 

HD 58343 

1.2 

± 

0.4 

2.6 

± 

0.7 

9.69 

± 

0.06 

1.6 

± 

0.4 

2.5 

± 

0.2 

11.84 

± 

0.02 

0.47 

± 

0.05 


3337 

± 

322 

HD 76341 

5.8 

± 

0.1 

3.1 

± 

3.1 

9.77 

± 

0.08 

6.3 

± 

0.8 

4.8 

± 

0.1 

12.12 

± 

0.02 

0.29 

± 

0.04 


2106 

± 

312 

HD 112272 

-17.7 

± 

0.4 

3.5 

± 

0.7 

9.88 

± 

0.08 

-16.9 

± 

0.1 

3.5 

± 

0.1 

12.16 

± 

0.02 

0.35 

± 

0.05 


2474 

± 

367 

HD 115842 

-14.8 

± 

0.5 

0.4 

± 

0.4 

9.7 

± 

0.08 

-14.9 

± 

0.2 

3.3 

± 

0.3 

12.24 

± 

0.02 

0.19 

± 

0.03 


1360 

± 

201 

HD 136239 

-52.6 

± 

0.2 

2.3 

± 

2.3 

9.93 

± 

0.08 

-52.3 

± 

0.9 

5.5 

± 

0.5 

12.60 

± 

0.03 

0.14 

± 

0.02 


1008 

± 

123 

HD 142758 

-43.8 

± 

0.1 

1.1 

± 

0.0 

9.86 

± 

0.08 

-44.2 

± 

0.1 

5.0 

± 

0.1 

12.22 

± 

0.02 

0.29 

± 

0.04 


2058 

± 

305 


-11.7 

± 

0.2 

1.4 

± 

1.4 

9.54 

± 

0.08 

-12.2 

± 

0.2 

3.0 

± 

0.2 

12.12 

± 

0.02 

0.17 

± 

0.03 


1240 

± 

184 


2.3 

± 

0.5 

2.8 

± 

0.6 

9.43 

± 

0.08 

2.4 

± 

0.1 

5.0 

± 

0.1 

12.01 

± 

0.02 

0.17 

± 

0.03 


1240 

± 

184 

HD 148937 

3.0 

± 

0.2 

1.7 

± 

0.3 

9.67 

± 

0.08 

2.9 

± 

0.5 

3.1 

± 

0.2 

12.26 

± 

0.02 

0.17 

± 

0.03 


1212 

± 

180 

HD 148379 

-18.1 

± 

0.5 

5.8 

± 

1.2 

10.35 

± 

0.08 

-18.5 

± 

1.4 

8.0 

± 

0.3 

12.50 

± 

0.02 

0.47 

± 

0.07 


3337 

± 

494 


-12.1 

± 

0.1 

0.4 

± 

0.4 

9.96 

± 

0.08 

-11.4 

± 

0.9 

5.5 

± 

0.7 

11.75 

± 

0.04 

1.07 

± 

0.10 


7646 

± 

738 

HD 148688 

2.8 

± 

2.0 

3.1 

± 

0.6 

9.98 

± 

0.08 

3.7 

± 

0.9 

5.2 

± 

0.2 

12.59 

± 

0.01 

0.16 

± 

0.03 


1157 

± 

202 

HD 152003 

2.3 

± 

0.2 

4.5 

± 

0.9 

10.26 

± 

0.08 

2 

± 

0.2 

3.0 

± 

0.7 

12.17 

± 

0.04 

0.81 

± 

0.08 


5800 

± 

560 

HD 154811 

-19 

± 

0.4 

0.8 

± 

0.8 

9.21 

± 

0.30 

-19.4 

± 

0.1 

3.4 

± 

0.5 

10.89 

± 

0.04 

1.38 

± 

1.13 


9850 

± 

8074 


6.7 

± 

0.2 

3.1 

± 

0.6 

9.96 

± 

0.08 

6.5 

± 

0.6 

3.6 

± 

0.1 

12.16 

± 

0.02 

0.42 

± 

0.06 


2975 

± 

440 

HD 154873 

-23.7 

± 

0.1 

2.8 

± 

2.8 

9.12 ± 0.3 

-23.7 

± 

0.1 

5.6 

± 

0.2 

11.64 

± 

0.02 

0.20 

± 

0.18 


1424 

± 

1289 


-12.8 

± 

0.1 

3.8 

± 

0.8 

9.95 ± 0.08 

-13 

± 

0.8 

4.6 

± 

0.1 

12.18 

± 

0.01 

0.39 

± 

0.07 


2776 

± 

486 


0.5 

± 

0.1 

3.8 

± 

0.8 

9.95 ± 0.08 

0.3 

± 

0.5 

3.1 

± 

0.1 

12.03 

± 

0.01 

0.55 

± 

0.10 


3921 

± 

686 


7.3 

± 

0.1 

1.7 

± 

0.0 

10.07 ± 0.08 

7.1 

± 

0.5 

2.8 

± 

0.1 

12.19 

± 

0.02 

0.50 

± 

0.07 


3576 

± 

530 

HD 155806 

-36.5 

± 

0.4 

2.0 

± 

2.0 

9.47 ± 0.30 

-37 

± 

0.1 

4.9 

± 

0.5 

11.47 

± 

0.04 

0.66 

± 

0.54 


4714 

± 

3864 


-2.4 

± 

0.5 

1.8 

± 

1.8 

9.47 ± 0.30 

-2.2 

± 

0.4 

2.2 

± 

0.8 

11.52 

± 

0.06 

0.59 

± 

0.43 


4202 

± 

3100 

HD 163800 

1.2 

± 

0.0 

6.5 

± 

6.5 

9.48 ± 0.30 

0.6 

± 

0.4 

5.0 

± 

0.1 

11.55 

± 

0.02 

0.56 

± 

0.51 


4013 

± 

3633 


7.0 

± 

0.0 

2.1 

± 

0.4 

10.39 ± 0.08 

6.3 

± 

0.5 

2.8 

± 

0.5 

12.44 

± 

0.02 

0.59 

± 

0.09 


4202 

± 

622 


t - b value fixed 

* - Sightlines with varying components in either Ca or Na 

For HD 29138, the (1) and (2) indicate the epoch in which these column densities were measure, ie, (1) was in 2002 and (2) in 2008. As 
none of the calcium lines varied in the other sightlines between epochs, these values apply for both 2002 and 2008 data. 
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Cold and dense sheets or filaments aligned along the line 
of sight could explain the column density variations. |Heiles| 
also suggested that observations of high proper motion stars 
should show variation in almost all sightlines, which is not 
the case in this study. Also, evidence for the low tempera¬ 
tures (~ 10 K) required by the Heiles model have not been 
found in the cold ISM, although Jenkins & Tripp (20111 
have detected a small but signihcant fraction of the diffuse 
ISM which is at high pressure. 

Studies of trace and dominant ions in the UV have in¬ 
stead indicated that the variations may be due to local dif¬ 
ferences in ionization rather than density (Danks et al.|200f 


Lauroesch fc Meyer||2003 Welty|[2007 1. However, this does 


not easily explain directly-observed variations in Hi (Frail 


et al. 1994 Diamond et al. 1989 Brogan et al. 20051, nor 


the existence of au-scale structure in maps of dust emission 
(Miville-Deschenes et al. 20101. Some authors have proposed 
alternative models, in which the observed variations are due 
to statistical fluctuations caused by a power spectrum of 
cloud sizes ( Deshpande|2000 20071. 

To date, only about 25 examples of temporal variation 
are known ( Lauroesch|2007 1. Two notable searches for vari¬ 
ation were those by [Danks et al.| ( |2001[ ) who studied ten 
stars behind the Vela supernova remnant and detected four 
sightlines which varied, and Smoker et al. (20111 who found 
evidence of variation towards three out of 46 stars spread 
across the southern sky - similar to the numbers of varying 
sightlines found in this study. 


4.1 Cloud electron densities and estimates of 
number densities and sizes 

At the low densities in the ISM, ionization states are deter¬ 
mined by the balance between photoionization and recombi¬ 
nation. For neutral and first ionization states in equilibrium 
this can be expressed as; 

r{Xi)n{Xi) = ar{Xi, T)n{Xii)n, (2) 

where F is the photoionization rate, is the radia¬ 
tive recombination rate coefficient, is the electron num¬ 
ber density and n{X) is the number density of species X 
( Welty, Hobbs fc Morton||2003 K 

Column densities of Cal and Call can be employed to 
estimate the electron density of the absorbing medium if it is 
in Ca ionization equilibrium and at a temperature of about 
100 K. If we assume a constant density along the line of sight, 
n{X) may be replaced with N{X), i.e. number density with 
column density, allowing us to rewrite equation (2) for the 
specific case of Cal and Call which becomes: 


FjCai) NjCai) 

Oir{Caj^T) N[Cajj) 


( 3 ) 


It is necessary to assume that the Cal and Call lines 
sample the same material. In the sightlines where the 6-value 
for Cal is smaller than that of Call, it is more likely that 
the former lies in denser regions, due to enhanced recombi¬ 
nation. Hence, in the following discussion, it is important to 
bear in mind the caveat that Cal and Call may trace dif¬ 
ferent regions (See Crawford|2002 for further discussion). 

Using the methods of Welty, Hobbs & Morton (20031 


and Smith et al. (20131 we have estimated the electron den¬ 
sity of components in those sightlines where we see varia¬ 
tion, and also for a sample without (where Cal and Call 
were present). These values are listed in table The sam¬ 
ple was chosen using the criterion that both Ca ll and Ca 
I components found in Hunter et al. (20061 and Smoker & 
Ledoux (in preparation) had comparable velocities (within 
Ikms^'^of each other). As these values can only be calcu¬ 
lated when Cal is visible, Tablej^is therefore heavily biased 
towards the densest material, where Ca I is most abundant 
and thus easiest to detect. 

Assuming that the electrons released through photoion¬ 
ization of atomic carbon dominate the electron density 
( Crawford||2002 \, we can use this to estimate a lower limit 
to the total number density of the sightline. Adopting a car¬ 
bon depletion of 6 0% implies a gas-p hase carbon abundance 
of 1.4 X 10~*nH ( Sofia et al.||l997 K Values of uh for each 
sightline where we have calculated rie are shown in Table 
These density estimates (at least those with nu ~> 10^ 
cm“®) are consistent with the filament structure predicted 
by |Heiles| ( [l9^ , as found in |Crawfor'd| ( |2002[ ). Typical va - 
ues of Uh are 10 - 100, 500, and 500-5000 cm for the diffuse 
atomic, diffuse molecular and translucent ISM respectively 
( Snow fc McCall||2006 1 . As can be seen from Tableall of 
our Uh values lie in or beyond the upper range. Typical val¬ 
ues stated in Welty (20071 of rie for the diffuse neutral gas 
are « 0.05-0.20 cm“'’, calculated from C i, Na i and/or K i. 
Again all but five of the densities found in Table are larger 
than this ’’normal” range. Although the values of rie and uh 
are high for the components in which we detect variations, it 
is interesting to note that they are in agreement with those 
for components that do not show variation. 

It has generally been found that electron densities es¬ 
timated from calcium tend to be higher than those derived 
from other species observed in the same sightline (e.g.| Welty 


et al. 1999 Sonnentrucker et al. 20031. Welty, Hobbs & Mor¬ 


ton (2003 \ have discussed this effect in detail. The cause of 
this effect is unknown ( [Welty, Hobbs fc Morton||2003[ ), but 
it has only been reported towards a handful of stars. In the 
absence of other tracers we cannot exclude this possibility. 
It is also possible that the derived densities are underesti¬ 
mated due to the aforementioned caveat that Cal and Call 
may trace different regions. If Call is mostly from an outer 
envelope, and Cal from a more dense region as discussed in 


Crawford (20021, then the Cal/Call ratio (and hence the 
density) will be greater than the value listed in Table 

We were able to calculate densities for 32 velocity com¬ 
ponents towards 19 stars, of which two components were 
found to include time-variable species (Sect. 3). High num¬ 
ber densities of > 10® cm“® are derived for all of the com¬ 
ponents. It appears that at least some regions of the diffuse 
ISM are at densities which are much higher than the typical 
values of 10® —10^ cm“®. If the assumptions leading to these 
values are correct, they lend cr edence to the hig h densities 
of (1 — 2) X 10® cm“® derived by Crawford (20021 and Smith 


et al (2013) for k Vel. It is also consistent with the detection 
of a small but significant population of diffuse ISM compo¬ 
nents which are at much higher pressure than the surround 
medium ( Jenkins fc Trip^|2011 |. Conversely, density esti¬ 
mates from C I by Jenkins & Tripp (20111 are usually sig- 
nihcantly lower than those estimated from Cal/Call (Welty 
et al. in preparation). 
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4.2 Biased Neutral formation in the ISM 


As in previous ISM studies [see] [for a review]citeplau07, 
changes in Nal are more often observed than in Call. [Lau^ 
roesch & Meyer (20031 estimate the minimum apparent 


length for the size of fluctuations observable in different 
species, being set by the distance an atom can travel be¬ 
fore It is ionised, i.e. the velocity divided by the photoion¬ 
ization rate. Using reasonable values for the velocity (es¬ 


timated from the ^-values) and photoionization rate, Lau 


roesch & Meyer (20031 estimated that Nal has an apparent 


scale of around 250 au, five times smaller than predicted 
for Call, hence explaining why typically tiny-scale struc¬ 
ture is more often seen in Nal. However, due to the higher 
photoionization rate of Cal, this species should show larger 
time-variability than Nal. On the other hand, the relative 
weakness of the Cal line means this is hard to test in prac¬ 
tice, although we have attempted to see if such variation 
would be detected in the Cal spectra where there is corre¬ 
sponding time-variability in Nal. 

In the current sample there are five sightlines for which 
we have multi-epoch data for Cal. Of these, one shows vari¬ 
ability in both Cal and Call (HD 29138) and three either 
have no observed Cal (HD 50896 and HD 112842), or have 
S/N insufficiently high to have detected the same percent¬ 
age change in column density as seen in Nal (HD 113904). 
The remaining target where time variability is seen in Nal 
(HD 163758) shows no Cal variation, and hence provides 
evidence against the biased neutral formation model. How¬ 
ever, we caution that this star is WR, so there is always 
the possibility that we are observing changes in the circum- 
stellar environment. Our results are therefore inconclusive 
in proving or disproving the biased neutral formation in the 


interstellar gas model of Lauroesch & Meyer (20031. 


5 SUMMARY 


We have observed 104 early-type stars at twin epochs to 
search for time-variable structures in the ISM over peri¬ 
ods of 5 to 20 yr. Our data set included 83 spectra with 
R ~ 80,000 from UVES and 21 spectra with R ~ 140,000 
from the 2.7m telescope at McDonald Observatory. Exclud¬ 
ing changes in absorption likely caused by eclipsing binaries 
or circumstellar discs, we find 5 cases within the UVES data 
of possible variation in the ISM. We have obtained follow¬ 
up third-epoch UVES and McDonald observations for the 
few objects which show possible change in their interstellar 
profiles. In only one case do we detect further variation. We 
expect this lack of variation, compared to [Watson fc Meyer| 
119961 who found differences in Na l towards the individ¬ 


ual members of all 17 of the binary or multiple star systems 
they observed (with separations 480-29000 au), is due to the 
much smaller scales probed in our investigation, i.e. from ~ 
0 - 360 au. Hence, although fine-scale structure seems ubiq¬ 
uitous on scale > 480 au it does not seem to be the case 
on scales < 360 au. We were able to calculate approximate 
densities, derived from the ratio of Cal and Call, in 32 
velocity components towards 19 stars, of which two compo¬ 
nents were found to include time-variable species (section 3 
and 4). Density estimates are consistent with the hlament 


structure predicted by Heiles 119971, as more recently found 


Crawford (20021. The majority of electron densities are 

de¬ 


higher than the typical values of 0.05 - 0.20 cm 


rived from other trace-dominant ratios (Welty 20071. The 
components in which we have found variation show no sig- 
nihcant difference in density calculations to those where we 
found no variation. Future work will include EEROS and 
HARPS data taken from the ESO archive to study an even 
larger sample of objects and investigate tiny-scale structure 
in terms of the power spectrum of the ISM. 
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APPENDIX A: POSSIBLE VARIATION IN 
TWO MCDONALD STARS 

Al HD 149757 Nal D variations 


HD 149757 Oph) is a runaway 09V B(e) star at a dis¬ 
tance of ~200 pc. The interstellar absorption-lines towards 
this star are very well studied, and the sightline often serves 
as an exemplar for modelling the ISM. Any variations in 
interstellar absorption towards this star are therefore po¬ 
tentially very important. Its displacement between the two 
epochs is 115.5 au. 

The Nal D absorption component at —5.6kms“^ shows 
a possible change between the two epochs separated by 20 
yr in Di. This variation is exaggerated when the McDon¬ 
ald data are compared with the models. However, when the 
original data (only available for Na D 5895.92 A, see Fig. 


Al ) were smoothed to allow the resolution to match that of 


the epoch 2 McDonald observations, the variation was less 
obvious. This implies a significant amount of the variation 
detected was due to the difference in resolution. Nonetheless 
this variation is still seen (however weak) in the compari¬ 
son with the epoch 1 data. This sightline has been studied 
in detail by Barlow et al. ( 1995| ) and so we have also pro¬ 
duced a fit using the velocity and column density values 
taken from Barlow et al. (19951, but with the instrumental 
&-value from the McDonald observations replacing those of 
UHRF. The resulting fit also indicates that the component 
at —5.6kms“^shows some possible variation. 

To investigate this variation further, we searched the 
ESO archives and found four observations with HARPS of 
Nal D towards this star, taken between 2006 and 2014. No 
further variations were detected, but from 2006 onwards the 
—5.6kms“^ absorption is saturated and so continued in¬ 
creases in column density might not be visible. Future inves¬ 
tigation of this component will need to examine the weaker 
near-UV lines of Nal at 3302 A, which are not included in 
the McDonald data. 

HD 149757 lies close on the sky to the LDN 204 fil¬ 
amentary cloud complex and ionizes the Hll region Sh 2- 


27. Radio observations of (Le Bourlot, Gerin & Per- 


aul^ 1989[ [Tachihara et al(]|2000 


[Liszt, Pety fc Tachihara 


20091 indicate that the cloud is accelerated by the radiation 


from C, Oph, causing compression of the molecular cloud and 
dense cores. This is supported by the CO observations of 
[Tachihara et ah] ( |2012[ ), who find several cloudlets of scales 
~1000 to 10,000 au which they suggest are cold, dense and 
overlapping. These structures could explain the apparent 
small-scale variations. 


HD 149757 



Figure Al. Overlay of WHK spectra (solid red line) and McDon¬ 
ald data (dashed blue line line) for HD 149757 (^ Oph) showing 
variations at —5.6kms“^ in the line profile for the Nal D line at 
5895 A. The residuals between epochs are plotted below. It was 
necessary to degrade the resolution of the data used by WHK. 


Their epochs overlap those probed by our McDonald spec¬ 
tra, but used a different instrument. They found that the 
column density of Nal at n = 18kms“^ (Heliocentric ve¬ 
locity) decreased by « 20% in 8 yr, but saw no change in 
Call. We confirm the variations in Nal found by [Lauroesch[ 
[fc Meyer] with the same component showing a similar varia¬ 
tion. This is evident from (Fig. A21 but, as was shown for HD 
149757, the difference in resolution between the spectra has 
caused some problems in the analysis of this sightline. Due 
to the limited spectral range available in the WHK spectra 
there is not much continuum visible on either side of the 
interstellar features and so normalisation would be difficult 
in this case. As noted by Lauroesch fc Meyer (20031, the 
continuum may have been fitted slightly too low, near the 
varying component at 23 kms”'^, however adjusting for this 
would only yield even larger differences and so cannot be an 
explanation for this variation. The McDonald data are plot¬ 
ted in the figure along with the WHK spectrum, degraded 
in resolution to allow for a more accurate comparison. This 
variation is apparent in the residual plot and so this may 
be showing quite strong variation in the ISM over « 20 yr. 
However, all components of this line vary between epochs. 
The most redward component varies in the opposite direc¬ 
tion to the others (believed to be affected by differences in 
resolution ) which helps to confirm the variation, but makes 
quantifying it difficult. This variation is over a scale up to 
184 au, the transverse distance moved by the background 
star. 0 


A2 HD 91316 Nal D variations 

HD 91316 (p Leo) is a well-studied pulsating variable star at 
a distance of 1350 pc, and is a frequent target for interstel¬ 
lar absorption-line studies. Lauroesch & Meyer (20031 pre¬ 
viously discovered variations in the Nal D absorption-lines 
between 1995 and 2003 due to stellar motion of ~ 12 au. 


® We are unable to reconcile the proper motion measured by van 
Leeuwen (2007) with the transverse distance of 12 au in 8 yr 


reported by [Lauroesch fc Meyer[ 
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HD 91316 



Figure A2. Overlay of WHK (solid red line) and McDonald 
spectra (dashed blue line line) for HD 91316 (p Leo) showing 
variations at 23kms“^ in the Nai D line profile at 5895 A. The 
residuals between epochs are plotted below. It was necessary to 
degrade the resolution of the data used by WHK. 




